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We provide an account of the static and dynamic properties of hard-core bosons in a one- 
dimensional lattice subject to a multi-chromatic quasiperiodic potential for which the single-particle 
spectrum has mobility edges. We use the mapping from strongly interacting bosons to weakly inter- 
acting fermions, and provide exact numerical results for hard-core bosons in and out of equilibrium. 
In equilibrium, we find that the system behaves like a quasi-condensate (insulator) depending on 
whether the Fermi surface of the corresponding fermionic system lies in a spectral region where the 
single-particle states are delocalized (localized). We also study non-equilibrium expansion dynamics 
of initially trapped bosons, and demonstrate that the extent of partial localization is determined by 
the single-particle spectrum. 



I. INTRODUCTION 

The localization of quantum particles and waves in dis- 
ordered media has been a key theme common to several 
fields of physics since Anderson's prediction of the effect 
half a century ago [lj. Among the different setups lead- 
ing to localization effects, quasiperiodic potentials with 
incommensurate periods [2 J provide a particularly attrac- 
tive realization, as their mathematical simplicity allows 
for various exact results. A one-dimensional (ID) sys- 
tem subject to such a quasiperiodic potential (Aubry- 
Andre potential [2]) has a localization transition at a 
finite value of the potential strength, in contrast to a 
random potential which in ID causes localization at in- 
finitesimal strength. Recently, pioneering experiments 
have explored the physics of bosons [SHS] and light 18] 
in such potentials. 

The single-particle spectrum acquires additional struc- 
ture upon modification of the quasiperiodic superlattice 
potentials e.g., when there are two superlattice potentials 
with different wavelengths each incommensurate with the 
lattice, or by other modifications of the basic Aubry- 
Andre (AA) potential [9H2T]. In these cases, different 
energy regions of the single-particle spectrum might be- 
come localized at different strengths of the potential, so 
that at some strengths of the potential there are both 
localized and delocalized eigenstates. Such a structure is 
known as a mobility edge |9J. Mobility edges are well- 
known to exist in truly random potentials, but only in 
higher dimension. 

The interplay between disorder and interactions is a 
prominent theme in the study of strongly correlated sys- 
tems. With the realization of quasiperiodic potentials 
hosting bosonic atoms jBHSJ? a natural question is the be- 
havior of interacting bosons in quasiperiodic potentials. 
Refs. [22-38J have studied the Bose-Hubbard model in 
the AA potential for nonzero interactions. Some of these 
works [32-35, 37J have used the infinite-interaction or 
hard-core limit of the Bose-Hubbard model, where mul- 
tiple occupancies are disallowed. The hard-core boson 
(HCB) model can be mapped onto free fermions, allow- 



ing numerically exact calculations for relatively large sys- 
tem sizes even in the absence of translation symmetry. 
Ground state, finite-temperature, and non-equilibrium 
properties of hard-core bosons in the AA potential have 
been analyzed in some detail in these studies [32 -35, 37J. 

In this article, we address ground- state and non- 
equilibrium properties of hard-core bosons in multi- 
chromatic quasiperiodic potentials (extended Aubry- 
Andre models) where the single-particle spectrum dis- 
plays mobility edges. Our main result concerning equilib- 
rium properties is a connection between the behavior of 
the many-body system (insulating or quasi-condensate) 
and the location of the Fermi energy of the corresponding 
free-fermion system (which is the chemical potential for 
both the corresponding free-fermion system and for the 
HCB). The system acts like a quasi-condensate (or insu- 
lator) when the chemical potential is in an energy region 
where the single-particle states are extended (or local- 
ized), irrespective of the nature of single-particle states 
lying below the chemical potential. This phenomenon is 
particularly puzzling in the case where the chemical po- 
tential is in a localized region but there are filled extended 
states at lower energies. One would expect superfluid 
properties due to extended states being occupied. We 
demonstrate this remarkable property, that the system 
behavior is determined by the location of the chemical 
potential, through the study of various quantities (quasi- 
condensate fraction scaling, off-diagonal order, entangle- 
ment entropy, etc). 

We also study the dynamics after release of an initially 
trapped cloud of hard-core bosons in multi-chromatic 
potentials. Expansion of initially confined many-body 
systems in potentials displaying mobility edges (higher- 
dimensional random potential) has been explored re- 
cently in two experiments [39| [40 ] . providing motivation 
for theoretical non-equilibrium calculations such as ours. 

We will start in Section [IT] with a description of 
the model and some single-particle properties of multi- 
chromatic potentials. Section [TTT] presents various ground 
state properties, highlighting the connection to the posi- 
tion of the chemical potential. Section |IV| discusses the 
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FIG. 1. Inverse participation ratio of all single-particle eigen- 
states, for the Aubry- Andre (AA) model and for the extended 
models EAA-1, EAA-2 and EAA'. Here i is the eigenstate in- 
dex in order of increasing eigenenergies. Except for AA, at 
intermediate potential strengths Vi both localized and delo- 
calized eigenstates coexist. Mobility edges are observed as a 
function of energy as one crosses boundaries between localized 
and extended regions. 



non-equilibrium expansion dynamics. 



II. MODEL AND METHOD 

We consider a system of hard-core bosonic atoms 
on a chain with L sites and open boundary conditions, 
described by the Hamiltonian 



H = 



jJ2{ b n b ^+ b Ul h n) + E^)^n (1) 



where b n and b\ are bosonic creation and annihilation 

operators, b n ,b\, =5 n ^ n >. The hard-core constraint is 

imposed through b 2 n = b^ = 0. We use several quasi- 
periodic potentials, V(n), specified below. 

Such a hard-core bosonic system may be regarded as 
the limit of the Bose-Hubbard model where the on-site 
Bose-Hubbard interaction is much larger than all other 
energy scales in the problem. This limit has the advan- 
tage of lending itself to exact calculations, through map- 
ping to a fermionic system via the Jordan- Wigner (JW) 
transformation. Introducing the fermionic operators c n 

(with |c n ,c^| = 8 n , n >) the mapping 



S(0,n- l)c n 



(2) 



where S{n,n') = Yl™ =n (1 — 2nj) is the JW string oper- 
ator, reproduces the HCB commutation relations. After 
the JW mapping the fermionic Hamiltonian can be ob- 
tained from Eq. replacing b n — » c n . The bosonic 
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FIG. 2. Left: Single-particle spectrum as a function of 
the quasi-periodic potential amplitude for the EAA-1 model. 
The color code indicates IPR values. The localization- 
delocalization transition occurs where many eigenstates 
bunch together; this bunching occurs at different Vi for differ- 
ent parts of the spectrum. Upper right: some representative 
wave functions for V\ = 2. 5 J (mobility edge region) for the 
states indicated with colored dots in left panel. Both localized 
and extended eigenfunctions are seen. Lower right: absolute 
value of the localized wave functions is plotted in log-linear 
scale showing that localization is exponential. 



one-body density matrix 



Pnn 



(3) 



can be computed from the fermionic one, p^ n , = 
using the approach described, e.g., in Ref. [41 J. 

For uniform ID bosons at zero temperature, the oc- 
cupancy of the lowest single-particle state (with momen- 
tum k = 0) scales as rik=o ^ V^b- While Bose Einstein 
condensation (BEC) is not strictly present in ID, the 
off-diagonal elements of the density matrix still develop 
a algebraic decay p^ n , oc _ 1 1/2 as \n — n'\ —> oo in 



T 



the thermodynamic limit. In spatially inhomogeneous 
situations, the quantity analogous to rik=o is the largest 
eigenvalue Ao of the single-particle density matrix. The 
eigenvalues are referred to as occupation numbers of nat- 
ural orbitals [4TH43] , The natural orbitals are the corre- 
sponding eigenvectors: 



p B <$> n 



An*?, 



with Ao > Ai > .... Quasi-condensation is signaled by 
the behavior Ao oc ^/Nb. Since we are dealing with ex- 
plicitly non-uniform systems, we choose to use this lan- 
guage (rather than momentum occupancies) in order to 
describe the presence or absence of quasi-condensation. 

We consider different types of quasi-periodic poten- 
tials. The simplest is a single-frequency cosine periodic 
potential with an irrational wavevector relative to the 
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lattice spacing, generally known as the Aubry- Andre po- 
tential [2J: 



i/L 



V(n) = V\ cos (27rqix n ) ; with q\ = 



a/5-1 



[AA] 

Here x n = an with a the lattice constant. The single- 
particle eigenstates of the AA model are all localized for 
V\ > 2 and all extended for V\ < 2. In this paper we are 
interested in extended Aubry-Andre models where the 
single-particle spectra have more intricate structure, in 
particular mobility edges. We introduce two potentials 
with two frequencies: 



V(n) = V\ cos (27rgix n ) + V 2 cos (2i\q 2 x 
The first model has the parameters 
V5-1 -1 



[EAA-1,2] 



V2 = \V U qi 2 

and the other has parameters 
0.7 , 



, q 2 = 3 Ql [EAA-1] 



V 2 = |Vi, gi —a , g 2 2<7i . 



[EAA-2] 



We also consider another type of modification of the AA 
potential: 



V(n) = Vicos(27rgi< d ); 
qi = — - — a , a d 0.7. 



[EAA'] 



These models EAA-1,2 US] and EAA' [12] are known 
examples from a large class of extended Aubry-Andre 
models whose single-particle spectrum possess mobility 
edges. 

As a function of the ratio V\jJ these models present 
the same qualitative features: for small values of V\j J 
all the single particle eigenvectors are extended in real 
space and for large values of this ratio all the eigenstates 
are localized. For an intermediate range of V\j J both 
localized and delocalized eigenstates coexist. The sin- 
gle particle spectrum is organized in rather well defined 
regions corresponding to localized or extended eigenvec- 
tors separated by sharply defined mobility-edges. The 
generic behavior of the localized-extended transition is 
illustrated in Fig. [T] where the inverse participation ra- 
tio IPR (ip) = r ^^L 2 is plotted for the single-particle 

eigenstates of of H. Note that for the A A model there are 
no mobility edges in the spectrum, instead, the localized- 
extended transition occurs for all eigenstates simultane- 
ously for the same value of V±/J. This non-generic fea- 
ture is due to the self-duality of the model |2J. Fig. [2] 
shows the single-particle spectrum as a function of V\ / J 
and illustrates the difference between (exponentially) lo- 
calized and extended states. 

Fig. [l] illustrates that the models EAA-1, EAA-2 and 
EAA' present the same generic behavior. We therefore 
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FIG. 3. Lower panel: quasi-condensate fraction given 
by Xo/\^Nb as a function of the filling fraction v — Nb/L. 
(Vi = 2.5J, EAA-1.) Quasi-condensation is observed when 
the chemical potential fi (y) lies in a spectral region with ex- 
tended states (in this case for v around half-filling, as seen 
from the top panel) . Otherwise Ao / \/Nb vanishes in the ther- 
modynamic limit Nb —> 00. Upper panel: inverse participa- 
tion ratio against eigenstate index. 



present our results mainly for one of the models (EAA-1). 
The expectation is that the concepts and results emerging 
from this work are generically valid for any potential for 
which the single-particle spectrum has mobility edges of 
this generic type. 



III. GROUND-STATE PROPERTIES 

In this Section we will present results characterizing 
the ground state properties of our system, as a function 
of the filling fraction v = Nb/L. The energy of the last 
filled JW fermionic level is the chemical potential of the 
system, /i(^). Our major result is that the system be- 
haves like an insulator or a quasi-condensate, depending 
on whether fi (y) lies in an energy region of localized or 
extended single-particle states, which we will denote re- 
spectively by Tii and E e . 

In order to characterize the ground state we consider 
the behavior of the natural orbitals occupation, the char- 
acteristic decay of the off-diagonal density matrix com- 
ponents, and the entanglement entropy of a subsystem. 
For definiteness we display results for the EEA-1 model 
with Vi = 5/2 J, where the single-particle spectrum has 
a well-defined intermediate delocalized region (E e ) sep- 
arated by mobility edges from higher and lower energy 
regions of localized states (Ej). 

When all the single-particle eigenstates are localized 
(ID random potential or AA model with V\ > 2), 
^o/V^6 ~^ [33 ; 37J, in contrast to the quasi-condensate 
behavior Aq ~ \fNb- In addition, in the localized (insu- 
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FIG. 5. (Vi = 2.5J, EAA-1.) (a) Lengthscale £ of expo- 
nential decay of the off-diagonal elements of p B , obtained by 



fitting with an exponential form p (x) 



-1*1/* 



(see text). 



The error bars are given by the fit to the numerical data. For 
the v values where quasi-condensation occurs, the behavior is 
algebraic (Fig. [4| and no £ values are shown, (b) Power-law 
exponent for off-diagonal decay p B (x) oc \x\~ a in the region 
where quasi-condensation is present. 



FIG. 4. Decay of the off-diagonal elements of p B for three 
different filling fractions in log-linear (upper row) and log-log 
(middle row) scales. (Vi = 2. 5 J, EAA-1.) Data points are 
shown for three different system sizes L to indicate the degree 
of convergence. The straight lines are fits. For v — 0.25 
and 0.75 [p (y) G Ej, no quasi-condensation], the decay is 
seen to be exponential; p B (x) oc Ae"^'^. For v — 0.5 the 
decay is well approximated by power law p B (x) oc a \x\~ a 
with a ~ 0.59. Lowest row: the natural orbital occupations 
A n . For v — 0.25 and 0.75 the distribution is step-like. For 
v — 0.5 it diverges for n— »0 in the thermodynamic limit. This 
can be seen in the inset where the black line corresponds to a 
divergence of the form A n oc c(n/L)~ K with n c± 0.3. 



lating) case, the off-diagonal correlations decay exponen- 
tially p B n , oc e ~\ n ~ n '\^ [331 E?J, in contrast to the quasi- 
condensate behavior p B , ex : hrm - In the A A model 

' nn \n— n'\ 1 

the two behaviors are seen for V\ > 2 an V\ < 2 at any 
filling fraction [33]. We will show that when the single- 
particle spectrum has mobility edges, either behavior can 
appear, depending on the filling. 

Fig. [3] shows the rescaled lowest natural orbital occu- 
pation Ao / \ /r Nb as a function of the filling fraction. The 
finite size scaling shows that if p (y) G £/, Ao/v^6 
in the thermodynamic limit, while Xq/^/NI approaches a 
finite value for p, (y) G E e . 

The stark difference between the p, (y) G E^ and 
pi (y) G E e cases is also observed in the off-diagonal ele- 
ments of the averaged one-body density matrix p B (xj) — 
Z Si ( see Figs. 4 pk. Here the average is taken 

to smoothen oscillations introduced by the quasi-periodic 
potential. For p(v) G E^ the decay is exponential 
p B (x) (X e~\ x \/£ whether for p(v) G E e one observes 
an algebraic behavior p B (x) cx We found the 

exponent a to be in the range a G [0.5,0.7] for our 
system sizes of (9(10 3 ), for EAA-1, EAA-2, EAA' mod- 
els. This is higher than in the bare tight-binding lattice 
(Vi = 0) case, for which a = 1/2 |41| . It is currently not 



known whether a would converge to 1/2 in the thermo- 
dynamic limit in multi-chromatic potentials with fillings 
corresponding to p (y) G E e . 

The occupation of the natural orbitals, displayed in 
Fig{4] (lower panel), shows that for pi (y) G E^ the dis- 
tribution of the A n 's is step-like, contrasting with the 
smooth decay for the p (y) G E e case. Step-like features 
are usually assigned to fermionic distributions. Here the 
natural orbitals are localized [37J and the bosons behave 
much like in a fermionic system as no localized orbital 
can accommodate more than one HCB. 

Another striking difference between p (y) G E^ and 
p (y) G E e cases is observed in the behavior of A n at n—>0. 
For p(y) G E e one observes a divergence A n oc (n/L)~ K . 
For the example displayed in Figj4j k ~ 0.3. However this 
exponent does depends on v and on the specific model 
of disorder and can even be larger the value k = 1/2 
obtained for = 0. For p {y) G E^ no divergence is 
observed. 

Finally, we consider scaling of the GS entanglement 
entropy of a subsystem D as a function of the sub- 
system size Ljj. The entanglement entropy is defined 
as Sd = — Ttd [pd lnAo], wnere Pd = Tr^ [p] and 
p = |^o) (^o| are many-body density matrices and Trp 
(Tr^) denote the trace over the degrees of freedom in 
D (in the complement of D). For HCB the expression 
for the entanglement entropy simplifies [44J to Sd = 
— Vi In vi + (1 — Vi) ln(l — Vi) where the v^s are the 
eigenvalues of the two-body density matrix of the JW 
fermions restricted to the subsystem D, pf^ with z, j G D. 
For piy) G £/, Sd saturates as Ld —> 00. However for 
p(v) G E e it behaves as Sd — flnLo for large subsys- 
tem sizes Ld with a prefactor c c± 1. This represents the 
well-known logarithmic correction to the "area law" in 
gapless one-dimensional systems [44j|45], with the pref- 
actor given by the central charge c = 1 in this case. Once 
again, the bosonic system behaves like a gapless quasi- 
condensate or like a gapped insulator depending on the 
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FIG. 6. Ground state entanglement entropy Sd of a subsys- 
tem D, averaged over the central position of the subsystem, 
for various subsystem sizes Ld- (L = 1000.) For the v val- 
ues where quasi-condensation does not occur [p, (y) £ Ej], 
Sd saturates to a constant value with increasing Ld- In the 
\i (v) G E e region, *Sd grows logarithmically with Ld- A log- 
arithmic fit (inset) shows that in this region the coefficient c 
approaches the value c — 1 which is the exact result for the 
case where no quasi-periodic potential is present. 



location of the "Fermi energy" of the corresponding free- 
fermion system. 



IV. EXPANSION DYNAMICS 

In this section we address the non-equilibrium expan- 
sion dynamics of an initially trapped cloud of bosons after 
the trapping potential is turned off at time t = 0. The 
motivation for studying this type of dynamical protocol is 
that expansion dynamics is particularly sensitive to the 
localization properties of the underlying single-particle 
states: particles will mostly fly off or mostly remain in 
the initial region depending on whether the initial overlap 
is dominantly with extended or localized single-particle 
states. This type of physics has also been explored in two 
recent experiments, which use higher-dimensional disor- 
dered potentials having mobility edges [39| 140] . 

The initial state (t < 0) is the ground state of the 
Hamiltonian with an added harmonic trapping po- 
tential, i.e., 



H(t) = H + e(-t)J2w(n)blb n 



(4) 



no being a site near the center of the chain. We moni- 

r 21 1 / 2 

tor the cloud size a (t) = {xjit)) - (xi(t)y and the 

fraction of localized atoms N\ oc (t)/N b as a function of 
time. Here, N\ oc (t) is defined as the number of atoms 
remaining in the support S of the initial density distri- 
bution, S = {i : (hi(t = 0)) ^ 0}. 

In the presence of a trapping potential and for = 0, 
a single dimensionless parameter p = N^yJW/ J controls 
the behavior of the atomic density in the large- N b limit 
[4T] , This means that p% n H> g (np/N b ; p) /N b as N b — >• 

00, where g (y; p) is the normalized density distribution. 
For p ~ 2.6 — 2. 7a Mott insulator region builds up in the 
middle of the trap, where the density is pinned to unity. 

As in the last section we focus on the model EAA- 

1. We study the dynamics for different values of V\ and 
p. We separate the case p — ^ 00, for which a simple 
interpretation of the results can be given, from the finite 
p that leads to more complex dynamics. 



where 0(...) is the Heaviside step function and 

W(n) = W(x n -x no ) 2 , (5) 



A. Large p 



We start by addressing the large p limit where, in the 
initial density distribution, all the atoms are in the Mott 
phase. In other words, in the initial state a central block 
of the chain is completely filled while the rest of the chain 
is empty. FigQa) shows snapshots of the time evolution 
of the initial density profile (left panel) for different values 
ofVi. 

Figs.[7jb,c) display the time evolution of the cloud size 
and the fraction of localized particles. For V\ < 3.8J, 
a (t) ex vt (for t —> 00), which means that at least some 
part of the cloud spreads ballistically. The spreading 
velocity v is a decreasing function of V\. The threshold 
value Vi ~ 3. St is that above which there are no more 
extended eigenstates in the spectrum, as can be seen from 
FigfTJb). For V\ > 3. 8 J, a (t) saturates and the width of 
the cloud remains bounded for large times. 

From Fig. [7^c) one can see that, at long times, the 
fraction of particles that remain within the region S is 
given by the fraction of extended states in the spectrum 
(horizontal dotted lines). This may be understood by 
considering the evolution in the eigenbasis of the non- 
trapped Hamiltonian (W n = 0) with which the system 
evolves for t > 0. In the asymptotic long time limit 
the fraction of delocalized atoms is proportional to the 
overlap of the initial state with the extended eigenstates 
of the non-trapped Hamiltonian. If the number of ini- 
tially occupied sites is not too small, the initial state has 
approximately an equal overlap with all localized and ex- 
tended eigenstates. The fraction of particles localized in 
S should thus equal the fraction of extended eigenstates 
in the spectrum. 
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FIG. 7. Expansion dynamics of bosons initially clustered in the center of the lattice as shown in top left. (Nb = 100, L = 1000.) 
In the other (a) panels, the density profile at t = 500J -1 is shown for V\ — 1J, 2.5J, 4J (extended, mobility-edge, localized 
regimes), (b) Root-mean-square width of the expanding cloud. When there is a fraction of extended states and since they have 
a non-zero overlap with the initial state, the evolution is always ballistic for large times a (£) ocvt, nevertheless the pre-factor v 
gets smaller with increasing V\ as the fraction of localized eigenstates increases. When all the single particle states are localized 
the cloud width saturates for large times, (c) Time evolution of the number of bosons in the central region S (see text). The 
dashed lines give the fraction of delocalized eigenstates in the single-particle spectrum. 



B. Arbitrary p 

Next, we address the time evolution starting from the 
ground state of a harmonic trapping potential charac- 
terized by a finite value of p. Fig. [8] shows the asymp- 
totic long-time value of the fraction of atoms localized 
in S as a function of p for different values of V\ and for 
Nb = 25, 50, 100. The convergence of our results for in- 
creasing A/5 shows that the physics of the large- limit 
is well-represented in our numerical simulations and also 
that, for these values, effects of the location of the 
quasiperiodic potential is sufficiently averaged over. 

For p = (no trap) and V\ > 1.2 J all the atoms are 
localized as the ground state is obtained by filling the first 
Nb single-particle levels, which are all localized for the 
EAA-1 potential in this filling region. For small p (very 
shallow trap) one can still find N^ localized single-particle 
eigenstates within the trapping length £ — N^a/p. This 
explains why all bosons remain localized for small but 
nonzero values of p. 

However, as p increases, the number of localized states 
available within a region of size £ starts to be smaller than 
N\y and the overlap with delocalized states is then finite. 



For large p the results of the last section are recovered 
and the fraction of localized atoms equals the fraction of 
localized states in the spectrum. 



V. SUMMARY & DISCUSSION 

We have presented a study of a strongly interacting 
many-boson system in extended Aubry-Andre models 
for which the single-particle spectrum displays mobil- 
ity edges, i.e., in situations where both extended and 
localized states are present in regions (E e and £j) of 
the single-particle spectrum. The bosonic system was 
treated by mapping to free fermions, so that numeri- 
cally exact calculations are possible for relatively large 
systems. We have shown through non-equilibrium calcu- 
lations that expansion dynamics can be used as a probe 
of the sizes of S e and £/ regions in the single-particle 
spectrum. 

For ground-state properties, the most striking result is 
that the properties of the many-body system (insulating 
or quasi-condensate) is determined solely by the location 
of the chemical potential p{y\ which has the interpre- 
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FIG. 8. Main panel: fraction of bosons remaining localized 
for asymptotically long times after switching off the trapping 
potential, as a function of parameter p. (L = 2000, EAA- 
1, various Vi and iV&.) For large p the asymptotic value is 
given by the fraction of delocalized eigenstates in the spec- 
trum (horizontal dashed lines). For small p all the bosons 
remain localized after releasing the trap. The inset displays 
the IPR (as in Fig. [T]); the colored lines correspond to the 
values of Vi/J of the main panel. Upper panels: the ini- 
tial density profile is displayed for Vi/J = (blue) and for 
Vi/J= 1 (gray) , for different values of p. 

tation of being the Fermi energy of the corresponding 
free-fermion system. We have illustrated this by tak- 
ing a representative model where, in a range of potential 
strengths, increasing the filling fraction v takes the chem- 
ical potential from a \i (y) G £j region at low densities, 
through a region of \i (y) G S e region at intermediate z/, 
to a /i (y) G E/ region at high densities. For different 
fillings z/, we presented the occupations of natural or- 



bitals (A n , particularly the quasi-condensate occupancy 
Ao), the off-diagonal decay of the single-particle bosonic 
density matrix, and the entanglement entropy scaling. 
These observables all show characteristics of gapped in- 
sulators in fillings for which /jl (y) G £/, and show char- 
acteristics of gapless quasi-condensates at v values for 
which \i (y) G E e . 

At present we lack a simple explanation for this re- 
sult that the location (relative to regions of the single- 
particle spectrum) of the Fermi surface of the correspond- 
ing free-fermion system determines off-diagonal proper- 
ties of strongly interacting bosons. While properties of 
fermionic systems are often described by features near the 
Fermi surface, the properties we have shown are bosonic 
rather than fermionic. Also, it is rather unintuitive that, 
when \i (y) G Ej, no quasi-condensate properties are seen 
even if a sizable fraction of Jordan- Wigner fermions oc- 
cupy extended states. 

Ref. [46J has observed, for another case (free fermions 
in disordered 2D Chern insulators) where the single- 
particle states has different natures in different parts of 
the spectrum, that certain properties depend solely on 
the nature of the eigenstates at the location of the Fermi 
surface. This loosely similar observation in a very dif- 
ferent context suggests that the situation is generic for 
systems where the single-particle spectrum contains re- 
gions of different nature. 

These results open up several questions that deserve to 
be addressed in future research. Most prominently, one 
would like to know if this dependence on the location of 
the filling fraction extends far beyond the hard-core limit 
for bosons with strong but finite interactions (e.g., Bose- 
Hubbard chain with finite U). This would be particularly 
intriguing because such a system does not map exactly 
to free fermions, so that a simple picture of filling up the 
fermi sea does not hold. More speculatively, one might 
also wonder if bosons in higher dimensional lattices with 
mobility edges also have properties determined solely by 
the location of the chemical potential, since in higher 
dimensions there is no connection to a fermionic picture. 
Another interesting direction would be to study finite- 
temperature properties. Finite temperatures might allow 
the system to explore the S/( e ) region even when the 
chemical potential is in the E e m region, if the chemical 
potential is sufficiently close to a mobility edge separating 
the two regions. 
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